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SUMMARY 

Ion temperatures were obtained in the SUMMA mirror device by observ- 
ing the Doppler-broadened charge-exchange component of the 667.8 and 
587.6 nanometer He lines in He plasma and the Ha and lines in H 2 
plasma. The second moment of the line profiles was used as the parameter 
to determine ion temperature. Corrections for magnetic splitting, fine 
structure, monochromator slit function, and variation in charge-exchange 
cross section with energy are derived and included. Even for constant 
cross section, no magnetic splitting nor fine structure, and infinitely 
narrow slit function, the line profile is not Gaussian, because the exci- 
tation results from a change-exchange process. Comparison is made with 
temperatures from a neutral particle analyzer. Electron temperatures 
were measured by the line ratio method for the corona model and were 
always much less than the ion temperatures. Correlations of ion and 
electron temperatures with plasma parameters are given. 


INTRODUCTION 

The purpose of the present work was to measure ion and electron 
temperatures and ion drift in the superconducting magnetic mirror appa- 
ratus (SUMMA) at NASA Lewis Research Center. The heart of this^facility 
is a steady-state magnetic-mirror device with crossed C- and B fields 
and superconducting magnets with the capability of producing a maximum 
midplane field of 4.9 tesla. The crossed and B fields cause the 
free electrons and ions of the plasma to drift azimuthally with approxi- 
mately the same velocity, corresponding to an ion drift energy of as much 
as several kilovolts. In addition, it appears that turbulence occurs 
which results in ion heating (ref. 1). The ion temperatures in such de- 
vices are almost always many times the electron temperatures. The SUMMA 
facility is unique in that it has the highest design field strength for 
its bore size of any existing superconducting magnet device. 

Ion and electron temperatures and ion drift were measured by emis- 
sion spectroscopy in the optical region from 412.1 to 728.1 nanometer 
wavelength. The ion temperatures were determined from the width of the 
charge-exchange neutral components of atomic spectral lines. The methods 
used in this spectroscopic determination of ion temperature and the com- 
parison of these temperatures with those determined in charge -exchange 
neutral detection is the principal subject of this report. An indication 
of the ion drift was obtained from the shift of the charge -exchange neu- 
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tral component of a spectral line relative to the narrow component asso- 
ciated with low energy neutral background particles • The electron temper- 
ature was found by taking ratios of intensities of spectral lines of 
helium and assuming a corona plasma model* 

Little previous work has been done on charge -exchange neutral Doppler 
broadening, which was first noted in aurora borealis* Recently Alexeff, 
et al. (ref. 2) observed the phenomenon in deuterium Balmer-a, -g, 
and -Y spectra from the turbulently heated plasma of the mirror device 
Burnout V and observed that . the most difficult problem here is 

that of appropriate unfolding of the line profile” (presumably in the 
presence of the narrow component and the slit function of the monochrom- 
ator). Hess (ref. 3) used charge -exchange neutral spectra to find the 
velocity distribution of a hot plasma beam from an r*f. -heated high- 
density plasma source. Sigman and Reinmann (ref. 4) observed a charge- 
exchange component of Hjg in a crossed-field mirror device. 

Much previous work has been done on determining free electron tem- 
perature from He line intensity ratios assuming a corona model* A review 
has been given by Sovie (ref. 5). Some high density effects have been 
discussed by Sovie in a later paper (ref. 6). The most successful appli- 
cations of the method have been made by Latimer, Mills, and Day (ref. 7) 
and Roth, Richardson, and Gerdin (refs. 8 and 9). 

Some previous spectroscopic measurements of ion and electron temper- 
ature and ion drift in SUMMA were made by Reinmann , Lauver, Patch, Posta, 
Snyder, and Englert (ref* 10) with a different electrode configuration. 

In the present work, the anodes and floating shields were changed to 
water-cooled copper, and the spectroscopic techniques were modified o The 
ion temperatures were obtained from stronger lines with no anomalies and 
less interference and were corrected for magnetic splitting, fine structure 
(where present) , monochromator slit function, and variation of charge- 
exchange cross section with energy. The ion drifts were obtained by taking 
averages of wavelength sweeps in both directions. The relative intensity 
calibration used in finding electron temperatures was improved. 

The scope of the present work was, nevertheless, somewhat restricted. 
No Abel inversions are included because the small viewing ports did not 
allow vertical scanning of enough of the plasma* Nor was an inversion 
method available for charge -exchange neutral spectra at specific wave- 
lengths. Transit time and high-density effects were not included. These 
and other effects are the subject of continuing research. 

Results are reported for pure helium 4 and pure hydrogen plasmas . 

Also, a mixture of 90 percent by volume H 2 and 10 percent He was used to 
measure electron temperatures from He line ratios. It is presumed that 
these measurements approximate the electron temperature in pure hydrogen 
plasmas * 

The authors are greatly indebted to Dr. Norman Tolk of Bell Tele- 
phone Laboratories for measuring the ratio of excitation cross sections 
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of the He 587,6 and He 667.8 nanometer lines. They are also indebted to 
Dr. Gordon Hammond of the U.S. Naval Ordnance Laboratory and Prof, James 
E. Bayfield of Yale University for helpful discussions. 


APPARATUS 
Magnet Facility 

The SUMMA facility (ref. 11) Includes four solenoidal magnets of 
vAiich the Inboard two were used for the present work (see fig. 1), The 
details of the magnet system have already been reported elsewhere 
(refs. 10, 12, and 13). The two inboard magnets alone produce maximum 
fields of 5.2 and 3,5 tesla at the mirrors and midplane, respectively. 


Plasma Test Section 


A schematic view of the plasma test section is shown in figure 1, 

The discharge chamber is a cylinder 3,75 meters in length and 36.6 centi- 
meters in diameter made from 304 stainless steel. A 10-inch diffusion 
pump at each end of the test section pumps the system to a base pressure 
of 5x10"^ torr. The large center port on one side was used for a neutral 


particle analyzer \diile the opposite port was used for spectroscopic ob- 
servations. 


Electrode Assembly 

The electrode arrangement is also shown schematically in figure 1. 
Jtore detail is shown in figure 2. Gas is introduced through the uncooled 
tungsten 12,7 millimeter-diameter cathodes to produce a hollow-cathode 
discharge. The water-cooled copper anode rings have 29 and 51 millimeter 
inside and outside diameters, respectively. They are electrically 
grovmded to the test section walls- The electrically floating shield pro- 
tects a cathode holder from heavy particle bombardment and helps inhibit 
arcing from the cathode holder to the anodes. The shield was constructed 
of copper and cooled by water. The anodes and cathodes are connected in 
parallel to a single d.c. power supply. 


Spectroscopic Diagnostic Equipment 

The apparatus for emission spectroscopy is shown in figure 3 . It 
was located at the magnetic midplane with the line-o£-sight always in the 
midplane. The vertical centerplane of the test section was focused on 
the entrance slit of the monochromator by means of a lens. Between the 
lens and the entrance slit, the beam was rotated 90° by means of a beam 
rotator consisting of totally reflecting prisms and mirrors. A stop was 
provided to reduce the vertical height of the beam. With this arrange- 
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ment, the area of the vertical centerplane viewed by the spectrometer was 
determined by the size and shape of the entrance slit and was approxi-* 
mately 2 millimeters high by 26 millimeters wide. Most of the light from 
the plasma came from inside a 54 millimeter plasma diameter centered on 
the magnetic axis* The vertical resolution of the optics was about 
4 millimeters after allowance for beam divergence over this distance. 

The monochromator, beam rotator, stop, and lens were rigidly secured 
to a table (fig* 3). The table was supported by a pivot and adjustable 
screw. With this arrangement, chords of the plasma at positions above 
and below the horizontal centerplane of the test section could be viewed. 
The distance of the chord from the centerline of the test section is 
designated y (fig* 3), 

The grating monochromator employed an f/8.6 ^ meter Ebert mounting 
with curved slits* The reciprocal linear dispersion was 1*6 nanometers 
per millimeter in first order* Detection was by means of a magnetically 
shielded photomultiplier with S 20 photocathode* The photomultiplier 
was cooled to about -30^ C by means of cold air from a vortex tube* 


Other Diagnostic Equipment 

A neutral particle analyzer (ref* 10) was located on the magnetic 
midplane on the opposite side from the monochromator and was used to 
measure energy distribution and mass of charge -exchange neutrals from 
the plasma* 

A floating electrical probe could be substituted for the neutral 
particle analyzer and was equipped for rapid traverse of the plasma. 


THEORY OF CHARGE-EXCHANGE DOPPLER BROADENING 
Line Shapes 

In most charge -exchange collisions involving electron exchange, the 
velocity of the projectile is little affected by the collision (ref* 14)* 
Consequently, for a charge -exchange process of the type 

Af + B A* + B+ (I) 

where the projectile is a positive atomic ion, the target B is a 

neutral atom or molecule. A* is an electronically excited charge -exchange 
neutral atom with the same nucleus as A"^, is B with one electron 
removed, and f indicates high velocity, a measurement of the A'*' tem- 
perature can be obtained by measuring the temperature of A* provided 
the A”^ temperature is much^greater^than the temperature of B. In low 
density plasmas in crossed £ and B fields this assumption is ful- 
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filled. The temperature Tj^ of A can be obtained by studying the 
Doppler broadening of the charge -exchange component of one of its spec- 
tral lines. 

The first step in Investigating charge-exchange Doppler broadening 
is to obtain the line shape function F, which is proportional to the 
emission intensity at a specific wavelength for an optically thin line 
and is normalized so that 



F(A)dA = 1 


( 1 ) 


where A is the displacement from the line center in terms of wavelength. 
It is assumed that has a Maxwell -Boltzmann velocity distribution 

(the loss cone of the mirror device and ion drift are neglected, and the 
cyclotron motion of the ions about magnetic field lines has small effect 
on the distribution at a given point in space) . It is shown in appendijc B 
that 



where o* is the charge-exchange cross section, E is ion energy in the 
laboratory system, and the other symbols are given in the LIST OF SYMBOLS 
(appendix A) , For the special case where the cross section o* is a con- 
stant, equation (2) becomes 


(3) 


If varies as v”^, equation (2) gives a Gaussian line shape function. 



A comparison of the well known Gaussian-Doppler line shape (ref. 15) 
and charge -exchange shapes from equations (2) and (3) is given in fig- 
ure 4 for He'*' at 3000 electron volts temperature. It can be seen that 
equations (2) and (3) give an appreciably wider line than the Gaussian- 
Doppler line shape. Also, the assumption of constant a* gives an appre- 
ciably different line shape than variable a*. 
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Second Moment Equations 

A closed form solution for may be obtained for the constant cr^ 

case if the Doppler broadening of the charge-exchange component of the 
spectral line is measured in terms of its second moment defined by 



(4) 


which is the average value of the square of the wavelength deviation of 
the line profile from its center. Use of the second moment also facili- 
tates corrections to given in appendixes C and D. 

It is shown in appendix B that if a* is constant 



(5) 


where the subscript 3 indicates that has been corrected for magnetic 

splitting, fine structure, and monochromator slit function; all of which 
would otherwise cause errors in Equation (5) may be solved for 

in closed form. 


It is also shown in appendix B that for the general case 




( 6 ) 


The corrections to <A y for magnetic splitting, fine structure, and 
monochromator slit function are derived in appendixes C and D. Let us 
assume we obtain the apparent shape of the charge-exchange component of 
the line by scanning the line with a monochromator and subtracting out 
the narrow component of the line (not due to charge exchange) and any con- 
tinuum or dark current. The apparent second moment can then be ob- 

tained with a digitizer attached to an electronic desk calculator (or, if 
necessary, with a digital computer). From equations (Cll) and (D15) 


n-fX-E 

i=l 


I.(X.. 
1 Oi 


Xq) 


w‘ 


8 In 2 


(7) 
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where the second term on the right Is the correction due to magnetic 
splitting and fine structure, and the last term on the right is the cor- 
rection due to the monochromator slit function. After the apparent sec- 
ond moment has been corrected by equation (7), equation (5) or (6) may be 
used to find (some such scheme as Interpolation of a table of 


and 


Ti 


has to be used with equation (6)). 


Selection of Spectral Lines 

Lines in the visible or near ultraviolet regions are the most con- 
venient experimentally. For hydrogen at 656.3 nanometer is the strong- 
est and has the least interference from H 2 lines. Hg at 486,1 nanometer 
is weaker and, under some conditions, has serious interference from as 
many as eight H 2 lines (ref, 10). For helium, the 587.6 and 667.8 
nanometer lines were found to exhibit the strongest charge-exchange com- 
ponents . 


Cross Sections 

The selection of appropriate charge-exchange cross sections for use 
in equation (6) is simplified by two facts: (1) Because the cross sec- 

tion appears in both the numerator and denominator, only the shape is im- 
portant - the absolute value i^ immaterial. (2) Below a certain projec- 
tile energy, when an ion collides with an atom, a collision complex is 
formed so that for identical nuclei the direct excitation and charge- 
exchange excitation cross sections are approximately equal. Hence, data 
on the sum of the two cross sections suffices. For He+, this occurs be- 
low about 15 kiloelectron volts (ref. 16), while for H"*" it occurs below 
about 5 kiloelectron volts (ref. 17) . 

In a hydrogen plasma, if the background gas is preominately H 2 , the 
following charge -exchange processes are important. 

+ H2(1so^) H^(n = 3) + H2(lsa) (II) 

hJ + H2(1so^) ^ H^(n = 4) + H^dsa) (III) 

Optical cross sections have been measured by Hess (ref. 18) for the H^ 
and Hg lines generated predominantly by processes (II) and (III) , respec- 
tively, for H'^ energies from 0.5 to 3 kiloelectron volts. These were ex- 
trapolated to lower energy by means of Massey's expression for the cross 
section of an inelastic atomic process (ref. 18). Hughes, Lin, and 
Hatfield (ref. 19) have measured these cross sections for H'*' energies 
from 5 to 130 kiloelectron volts. Hughes, et al. (ref. 20) have used ex- 
perimental 3s, 3p , and 3d cross sections to synthesize Improved H 
cross sections from 10 to 100 kiloelectron volts. It is Impossible to 
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fair Hess’s absolute cross sections into those of references 19 or 20 in 
a reasonable manner. Hesses work contains several possible experimental 
errors. Consequently, before fairing into cross sections from refer- 
ences 19 or 20, it was necessary to multiply all Hesses cross sections by 
a factor of 23, 

In a hydrogen plasma if the background gas is appreciably dissoci- 
ated, the following charge -exchange processes may be important, 

hJ + H(n = 1) ^ H^(n = 3) + H+ (IV) 

+ H(n = 1) H^(n = 4) + h'*' (V) 

No experimental or theoretical cross sections are available for these 
processes. For estimates, the total cross section for the process 

+ H(n = 1) ^ H^(n = 2) + (VI) 

was obtained by adding the 2s and 2p cross sections of Bayfield 
(ref. 21) and Young (ref. 22), respectively. The sum was extrapolated 
to smaller energies by Massey's expression for the cross section of an 
inelastic atomic process (refs. 18 and 23) with the internal energy dif- 
ference between the initial and final states AE ■ 10.20 electron volts 
and the effective interaction range a = 6.63xl0“8 centimeters. The 3s, 
3p, 3d, 4s, 4p , and 4d cross sections were then estimated from the 
n = 2 cross section by assuming the n and % dependence derived by 
May (ref. 24). The H optical cross section o was then calculated 
from (ref. 20) “ 

“c. ■ °3e % + “3d 

where the constant is a branching ratio. Similarly, the H„ optical 
cross section was calculated from ^ 

Op = 0.5842 + 0.1190 +0.7456 (9) 

where the three constants are branching ratios. 

For helium the most important charge-exchange processes are 

He^ + He(ls^ ^S) -»■ He^(3d ^D) + He'*’ (VII) 

He^.+ He(ls^ ^S) ^ He^(3d ^ + He'^ (VIII) 

for the He 587,6 and 667.8 nanometer lines, respectively. 
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The cross section of process (VII) has been measured by Muller and 
de Hear (ref- 25) between 1 and 150 kiloelectron volts. The total cross 
section for direct excitation and charge -exchange excitation has been 
measured by Dworetsky et al. (ref. 26) between 60 and 5000 electron volts 
on a relative scale. For this report the two sets of cross sections were 
normalized between 1 and 5 kiloelectron volts. 

The cross section of process (VIII) has been measured by Muller and 
de Heer (ref. 25) between 8 and 150 kiloelectron volts. The total cross 
section has been measured by Dworetsky et al. (ref. 26) between 60 and 
5000 electron volts on a relative scale. This leaves a region from 5 to 
8 kiloelectorn volts that has not been measured. More recent work by 
Tolk (ref, 27) disclosed that at 4 kiloelectron volts the total cross 
section for the He 667.8 nanometer line was 0.491 times the total cross 
section for the He 587.6 nanometer line. Hence, the Dworetsky, et al. 
(ref. 26) relative cross sections for the 667,8 nanometer line were mul- 
tiplied by 27x10“^® before fairing them into the charge -exchange cross 
sections of Muller and de Heer (ref. 25). 


Calculated Second Moments 

Second moments were calculated by equation (6) using cross sections 
from the previous section and lower and upper limits of integration of 
theoretical threshold and 1 megaelectron volt, respectively. The results 
are given in table I for the two best lines, assuming an H 2 backgrotmd 
gas for H(j. The results for in table I were based on the high-energy 
cross sections of Hughes, et al. (ref. 20), but did not change appreci- 
ably below a temperature of 1 kiloelectron volt (the region of Interest 
in this report), when the cross sections of Hughes, Lin, and Hatfield 
(ref. 19) were substituted. 

Second moment calculations for were also performed assuming an H 
background gas. As can be seen from figure 5, the assumption of an H or 
H 2 background gas had little effect on the calculated temperature for a 
given value of second moment. However, in view of the estimated cross 
sections for processes (IV) and (V) , this conclusion should be considered 
tentative . 


RESULTS AND DISCUSSION 
Test Conditions 

The principal plasma parameters were the midplane magnetic flux 
density, mirror ratio, total current for both sets of electrodes in par- 
allel, electrode voltage, and gas pressure. The midplane magnetic flux 
density was varied from 1.06 to 2.49 tesla. The mirror ratio varied from 
1.52 to 1.57. 
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Conditions for hydrogen . - Total current amounted to 0,60 to 
1.15 ampere. Electrode voltage was 7.0 to 20.4 kilovolts. Gas pressure 
was 6.7x10“^ to 1.12x10“^ torr; 

Conditions for helium . - Total current was 0.45 to 1.05 ampere. 
Electrode voltage was 4.9 to 19.9 kilovolts. Gas pressure was 1.04x10"^ 
to 2.02x10"^ torr. 

Conditions for mixture . - The mixture, which was used to obtain elec- 
tron temperatures in hydrogen, consisted of 90 percent H 2 and 10 percent 
He by volume. Total current was 0,84 to 0.86 anqjere. Electrode voltage 
was 11.9 to 18,8 kilovolts. Gas pressure was 1.13 to 1.24x10"^ torr. 


Ion Temperatures 

H ion temperature . - A typical recording of the line is shown in 
figure 6 for y = 0. The narrow central peak was caused by Franck -Condon 
neutrals and electron excitation of other neutrals and is of no interest. 
Since it had an approximately Gaussian line profile, it was easily faired 
out as shown. In some runs the H 2 657.2 nanometer line appeared (fig. 6). 
This was also faired out. The remaining curve was used to find the appar- 
ent second moment. The Hg line was similar. 

Equations (7) and a more extensive version of table I were used to 
find ion temperatures from the Hj;, and Hg lines. The sum of the second 
moment corrections for magnetic splitting, fine structure, and mono- 
chromator slit function never exceeded 11 percent. In general, the tem- 
peratures derived from and Hg, assuming H 2 background gas, were not 
exactly the same for y = 0 at identical experimental conditions, as fig- 
ure 7 shows. On the average, Hg gave temperatures that were 71 percent 
of those obtained from Ha. This is in contrast to the deuterium results 
of Alexeff et al. (ref, 2) who found that and Dg gave essentially the 
same ion temperatures. The cause of the discrepancy is not known, but 
calculations are in progress to see if the two-step process. 


e^ + H(n = 1) e + H(n = 4) 

(IX) 

H+ + H(n = 4) ^ H^(n = 4) + 

(X) 


which should be significant at high H'*' densities, might be responsible. 
Process (X) is a resonant charge-exchange process in contrast to proc- 
esses (II) to (V), 'Which are nonresonant. This affects the energy de- 
pendence of the cross section , as shown in figure 8. Moreover, the cross 
sections for processes such as (X) scale as n^ (ref. 28), whereas the 
cross sections for processes such as (II) to (V) scale as 1/n^ (ref, 17), 
Thus, it appears plausible that for Hg at sufficiently high ion densi- 
ties, process (X) might predominate over processes (III) and (V), but for 
the same conditions processes (II) and (IV) might predominate for H^j. 

This would make the charge-exchange component narrower than if proc- 
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esses (III) and (V) predominated (see fig. 8). 

Process (X) is not the only possible explanation for the discrepancy 
between Hqj and Hg temperatures. We have not considered the charge- 
exchange processes (refs. 2 and 18) 

H^^dsa) + H2(l30^) ^ Hj(n) + H(n') + hJ (XI) 

+ H^dsa) ^ Hj(n) + (XII) 

However, for process (XI), the cross sections for and Hg have about 
the same energy dependence (ref, 18) and scale as 1/n^ , where b is 
greater than 3, If the ion density were high enough, process (XII) might 
be important. Setting up to measure the electron density in SUMMA is in 
progress. 

Another possible cause of the discrepancy between and Hg is the 
longer radiative lifetimes of the n = 4 states of H. This would allow 
the n = 4 charge -exchange neutrals to travel further and, therefore, 
spread out more before radiating. The fast charge -exchange neutrals 
would spread out further than the slow ones, so Hg would appear narrower 
than for a homogeneous plasma. However, the effect shoxild be more pro- 
nounced at high ion temperatures^ which is not the case (see fig, 7). 

The experimental results are averages across a plasma diameter. It 
would be more convincing if they could be Abel inverted to give ion tem- 
perature as a function of radius. However, at a specific wavelength, the 
radiation at a point is not isotropic due to the Doppler effect, so it 
does not appear that an Abel inversion is applicable at individual wave- 
lengths , 

Two correlations were used to relate ion temperature to plasma 
parameters for y - 0, In figure 9, temperature derived from 
assuming H 2 background gas, is plotted against the ratio of power input 
to magnetic flux density. Less scatter is obtained if temperature is 
plotted against power input divided by pressure and magnetic flux density, 
as shown in figure 10 . 

A paired comparison of ten^eratures derived from and from a 
neutral particle analyzer is shown in figure 11, Each pair of temperatures 
were taken at the same time for the same plasma volume. The data are for 
y = 0, and the background gas was assumed to be H 2 . The agreement is 
satisfactory. 

ion temperature , - A typical recording of the He 587,6 nanometer 
line is shown in figure 12, The He 667,8 nanometer line was similar. 

The sum of the second moment corrections for magnetic splitting, fine 
structure, and monochromator slit function never exceeded 10 percent. 
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In general, the ion temperatures derived from the He 667,8 and 
587.6 nanometer lines were not exactly the same for y = 0 at identical 
experimental conditions, as figure 13 shows. On the average, the He 
667.8 nanometer line gave ion temperatures that were 65 percent of those 
from the 587,6 nanometer line. This may be due to the cross sections 
used. When He 587.6 nanometer charge-exchange cross sections were used 
for both lines, the temperatures agreed to within 3 percent. It will be 
recalled that the shape of the low-energy charge— exchange cross sections 
for both lines was obtained based on the assumption that the charge- 
exchange and direct excitation cross sections were equal. This may not 
be sufficiently accurate. Also, no cross sections for the He 667.8 
nanometer line were available between 5 and 8 kiloelectron volts and so 
were estimated. 

A more likely explanation for the lower ion temperatures from the 
He 667.8 nanometer line is obtained by considering the two-step processes 

rHe"^ + He ^ He^ + He'*’ (XIII) 

\^e^ + He^ e + He* (XIV) 


r e^ + He -V e + He* (XV) 

IjieJ + He* ^ He* + He'*' (XVI) 

Order of magnitude calculations indicate that these processes are some- 
^diat less important than processes (VII) and (VIII), and they cannot be 
included rigorously without knowledge of the free electron density and 
the geometry of the plasma . However , a conclusion may be drawn by 
examining the cross sections of processes (VII) , (VIII) , and (XIII) to 
(XVI). Process (VIII) has about half the cross section of process (VII) 
for typical energies, but processes (XIV) and (XV) have about the same 
cross sections for the He 587.6 and 667.8 nanometer lines. Process (XIII) 
is identical for both lines • There is no reason to believe that the 
cross sections for process (XVI) for the two lines are not essentially 
the same. Hence, processes (xni) to (XVI) should contribute relatively 
more to the 667.8 nanometer line than to the 587.6 nanometer line. Since 
processes (XIII) and (XVI) are resonant, they would make the 667.8 
nanometer line narrower (see fig. 8) and, if neglected, cause too low a 
temperature to be inferred. 

The discrepancy in ion temperatures cannot be explained as a radia- 
tive lifetime effect because the radiative lifetimes of the upper states 
of the two lines are essentially the same. 


Two correlations of ion temperature with plasma parameters were made. 
In figure lA, ion temperature for y = 0 is plotted against ratio of 
power input to magnetic flux density. Less scatter occurs in figure 15, 
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where ion temperature is plotted against power input divided by pressure 
and magnetic flux density. 

A paired comparison of ion temperatures from the He 587.6 nanometer 
line and from a neutral particle analyzer is made in figure 16. The 
cause of the discrepancy is unknown* 

Overview . - The measurement of ion temperature from the charge- 
exchange co^onen t of neutral atomic lines was found to be a simple pro- 
cedure allowing direct data reduction with a minimum of assumptions. It 
does not disturb the plasma. The temperatures so obtained are in reason- 
able agreement with temperatures measured with a neutral particle ana- 
lyzer. Accuracy can be improved by better models, more accurate cross 
sections, and Inclusion of additional char ge -exchange processes. 

A comparison of ion temperatures obtained rigorously (eqs^, (6) 
and (7)) with ion temperatures that would be obtained by erroneously 
applying the usual Gaussian Doppler half-width equation (eq. 4-117 of 
ref. 15) is interesting. If the charge -exchange process has a constant 
cross section, the erroneously obtained temperatijres will be 71 percent 
too high. If the cross section is a function of energy, the erroneously 
obtained temperatures will be roughly 50 to 240 percent too high for the 
plasmas of this reports 


Ion Drift 

If the guiding centers of the ions drift with a velocity component 
parallel to the line of sight, the entire charge-exchange component of a 
neutral atomic line will be shifted relative to the narrow central com- 
ponent because of the Doppler effect. This follows since the central 
component is due principally to electron collision excitation of neutral 
atoms, which are not expected to drift appreciably. Hence, the shift of 
the char ge -exchange component gives an indication of the drift velocity 
of the ions . 

Alexeff et al . (ref. 2) found that the observed shift also depended 
on the distances along the line of sight to the front and rear walls of 
their machine due to cascading and transit time effects. However, this 
spurious shift could be eliminated if the distances to the front and rear 
walls were equal . Shift measurements were made in SUMMA in hydrogen and 
helium, and no spurious shifts were observed despite unequal wall dis- 
tances. This was presumably due to the large inside diameter of the test 
section, which allowed decay before the neutrals struck the walls. 

Results of measurements of wide -component shifts as a function of 

y are shown in figure 17. These shifts are some sort of average across 
a chord, so it is difficult to interpret them quantitatively. However, 
if the shifts are interpreted as due to ^ x drifts, then they are 

approximately consistent with floating probe voltage measurements and cal- 
culated ^ fields. Both diagnostics indicate that the ^ field reverses 
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direction at roughly 3 centimeters. Figure 2 shows the probable reason 
for this. Between the anode and cathode B field lines » points inward. 
However, B field lines tend to allow a "short" between the cathode and 
the floating shield. The floating shield is then electrically connected 
by othe^ B field lines to the region outside the anode B field lines. 
Hence, C points outward from the anode B field lines. Thus, the plasma 
inside the anode B field lines rotates one direction and the plasma out- 
side the anode B field lines rotates the other direction. 

Similar shift results were obtained from the He 587.6 nanometer line. 


Electron Temperature 

The technique, used for these measurements has already been reported 
(ref, 10) , The method of relative intensity calibration was improved in 
the present work by using a larger tungsten strip lamp located on the 
centerline of the test section when in use and by taking into account, the 
wavelength dependence of the correction for reflection from the lamp 
envelope. 

To reduce the He line ratios to free electron temperatures, inten- 
sity ratios involving Maxwellian averages of optical cross sections are 
needed. The most reliable values for the 412.1, 443.8, 471.3, and 
504.8 nanometer lines have been given by Latimer, Mills, and Day (ref. 7), 
although one must note that their "intensity ratios" are really ratios of 
numbers of photons emitted. For the He 706.5 and 728.1 nanometer lines, 
cross sections have been given by Jobe and St. John (ref. 29) for low 
energies and by Moussa, de Heer , and Schutten (ref. 30) for high energies. 

All the following results are chordal averages for y = 0 (see 
fig. 3). 

Electron temperatures in He . - Ratios of He 412.1 and 443.8 nanometer 
line intensities were measured for a variety of conditions and used to de- 
termine the electron temperatures. The results are shown in figure 18 and 
are in. good agreement with reference 10. 

Ratios of He 504.7 and 471.3 nanometer line intensities were meas- 
ured for three conditions. They yielded electron temperatures an average 
of 60 percent higher than the He 412.1 and 443.8 nanometer lines. This 
compares with 29 percent higher in reference 10 and 27 percent higher on 
the centerline in BURNOUT VI after an Abel inversion (ref, 31) . 

Ratios of He 728.1 and 706.5 nanometer line intensities were measured 
for three conditions. They gave electron temperatures from 22 to 33 eV, 
an average of 13 percent higher than the He 443.8 and 412.1 nanometer 
lines . 

The disagreement of the temperatures from the various line ratios 
can be explained by applying Sovie's first criterion (table I of ref. 6) 
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to the singlet S depopulation rates using up-to-date cross sections 
(refs. 7, 29, and 30). This showed that there was appreciable depopula- 
tion of the n = 4 and n = 5 levels by electron collision deexcitation. 
Hence, the teinperatures from the n = 3 lines (728.1 and 706.5 nm) are 
the DQost reliable (ref. 6), 

Electron temperatures in a mixture . - It is impossible to obtain 
electron temperatures from H line ratios because the dependence of their 
excitation cross sections on electron energy is too nearly the same. Con- 
sequently, 10 percent He was added to the hydrogen. Only four suitable 
He lines were sufficiently strong: 471.3, 504.8, 706.5, and 728.1 

nanometers. Care was exercised to separate the He 471.3 and H 2 471.4 
nanometer lines and also the He 728.1 and H 2 728.0 nanometer lines. 

Electron temperatures were obtained from two ratios. Ratioing the 
728.1 and 706.5 nanometer line intensities gave temperatures from 13.0 to 
15.8 electron volts. These tenperatures would fall well below the temper- 
atures for pure helium if plotted in figure 18. This is what one would 
expect (ref. 32). Ratioing the 504.7 and 471.3 nanometer line intensi- 
ties gave temperatures about 2.5 times as high. 

The reason for the temperature discrepancy between the two pairs of 
lines is not the same as for pure helium. Applying Soviets first cri- 
terion gave no evidence of depopulation of the singlet 4s level by 
electron collisions. Instead, we believe the reason for the discrepancy 
is depopulation of the upper state of the 471.3 nanometer line by colli- 
sions of the second kind . 

He (4s ^S) + H^fClsa) ->■ He(2p ^P) + H(^S) + h'*' (XVII) 

The dissociation energy for the ground vibrational state of H 2 is only 
0.0208 electron volts more than the energy available when He drops from 
the 4s 3s to the 2p 3p state. Consequently, process (XVII) should 
have a large cross section (ref* 33). For deexcitation by process (XVII) 
to be just as Important as radiative deexcitation, process (XVII) would 
have to have a cross section of 4.9xl0~34 square centimeters for 10^3 
ions per cubic centimeter and an temperature of 1000 electron volts. 
This is a reasonable cross section if one considers the large orbital of 
the n = 4 state of He. A relatively large amount of at about this 
temperature has been observed in SUMMA with the neutral particle analyzer 
(ref. 10). The process analogous to process (XVII) but for singlet helium 
atoms would be expected to have a much smaller cross section because the 
energy deficit is 0.1951 electron volts instead of 0.0208 electron volts. 

In view of the above discussion, electron temperatures from the 
706.5 and 728.1 nanometer lines are more reliable. However, He 728.1 and 
H 2 728.0 nanometer lines are so close together that it might be more 
practical to use the He 706.5 and 504,8 nanometer lines to obtain elec- 
tron temperatures in the future. 
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COHCLUDING REMARKS 

Derivations are given which show that the determination of ion tem- 
perature from the charge-exchange component of a spectral line of a neu- 
tral atom can be put on a more quantitative basis if the second moment of 
the line profile is used as a measure of the ion temperature. Specifi- 
cally, use of the second moment makes it feasible to take magnetic split- 
ting, fine structure, monochromator slit function, and the energy depend- 
ence of the charge -exchange cross section into account. If one errone- 
ously applies the usual Gaussian Doppler half-width equation to a charge- 
exchange line component with a constant charge -exchange cross section, 
the calculated temperature will be 71 percent too high^ If the charge- 
exchange cross section varies, the similarly calculated temperature can 
easily be from 50 to 240 percent too high. 

The results were applied to measurements of ion temperature in SUMMA. 
It was concluded that the Balmer a line and the He 587,6 nanometer line 
were most suitable for this application in hydrogen and helium plasmas, 
respectively^ temperatures from 169 to 955 electron volts and He“^ 
temperatures from 666 to 5151 electron volts were obtained from these 
lines. In both plasmas,, the ion temperature increased with the parameter 
P/pB, where P is power input, p is gas pressure, and B is magnetic 
flux density. Ion temperatures so derived are averages across the plasma 
and are believed accurate to within 30 percent for and 35 percent for 
He"*", This is substantiated by a paired comparison with neutral particle 
analyzer results obtained at the same time for the same plasma volume « 

It should be possible to improve this accuracy by the incorporation of 
more accurate charge -exchange cross sections, a more sophisticated plasma 
model, and additional charge -exchange processes. 

Measurements of the shift of the charge -exchange component of a spec- 
^ral^line were also made. Results were approximately consistent with an 
C X drift velocity of the ions and indicated that the inner portion 

of the SUKMA plasma rotates one direction and the outer portion, the op- 
posite direction. For more quantitative results, a more sophisticated 
plasma model is needed o 

Electron temperatures were measured by the helium line ratio method 
assuml-Ug a corona model. Different line ratios gave somewhat different 
temperatures. The probable causes of this were diagnosed. It was con^ 
eluded that the ratio of the He 706.5 and 728.1 nanometer lines gave the 
most reliable temperatures in pure He and in H 2 ~He mixtures. This ratio 
gave temperatures between 22 and 33 electron volts for helium and between 
13 and 16 electron volts for a mixture of 90 percent hydrogen and 10 per- 
cent helium used to simulate pure hydrogen. 



17 


A 

a 

b 


c 

D 

E 

AE 



f 


k 

Z 

M 




P 

P 

S 

T, 

1 

->■ 

v,v 


APPENDIX A 
LIST OP SYMBOLS 

2 2 

constant in characteristic function, Me /2kTAQ 
effective interaction range 
magnetic flux density 

principal quantum number exponent for process (XI) cross section 
dependence 

velocity of li^t 

constant in characteristic function, ^ 
energy 

internal energy difference between initial and final states 
electric field 
line shape function 

line shape function observed at exit slit 
degenerate hypergeometric function 
Maxwell -Boltzmann velocity distribution function 
integrated intensity of component i 
Boltzmann * s constant 
azimuthal quantum number 
mass of ion projectile 
principal quanttim numbers 
power input. to plasma 
gas pressure 

slit function of monochromator 
ion temperature 
absolute velocity 


V M 
2kT 
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r-component of v 
z -component of v 


W 


full half width at half intensity of S in wavelength units 




difference between wavelength of light and monochromator 
setting, ^ 

second moment of S(x) 


y 

Y 

e 


distance from line of sight to centerline of test section 
(see fig. 3) 

2 2 2 

wavelength difference parameter. Me A /2 Aq 
polar angle 




displacement from the center of the charge -exchange compo- 
nent of a spectral line, X - 

difference between wavelength set on monochromator and 
center of line, Xj^ - Xq 


apparent second moment of charge-exchange component of 
spectral line 



corrected for magnetic splitting and fine structure 
corrected for monochromator slit function 


X 


wavelength of light 


m 


^Oi 


monochromator setting 

wavelength of center of charge-exchange component of spec- 
tral line 

wavelength of center of i^^ component due to splitting 

charge -exchange cross section for excitation to electroni- 
cally excited state 




Optical excitation cross sections 


^3s*^3p’^3d excitation cross sections for upper state given as subscript 
^4s ^^4p *^4d 


characteristic function for F 



characteristic function for S 


characteristic function for F 

e 

azimuthal angle 

Fourier transform argument 
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APPENDIX B 

DERIVATIONS OF LINE SHAPES AND SECOND MOMENTS FOR DOPPLER- 
BROADENED CHARGE-EXCHANGE NEUTRAL SPECTRA 
We consider the charge -exchange reaction 

A j + B A* + b'*’ 

and assume Aj has a Maxwell -Boltzmann velocity distribution and B is 
stationary. Davidson (ref, 34) has given the Maxwell -Boltzmann velocity 
distribution function f in polar coordinates. 

.3/2 .+Iv2/2kl^ 

f de d<j) dv = ^2ikf ) ^ ^ 

where 6 and 4> are the polar and azimuthal angles, respectively, be- 
tween V and the line of sight, v is absolute velocity, M is the mass 
of A"^, and k is Boltzmann's constant. The probability of Ax having 
its velocity vector in d9 d4> dv is directly proportional to f and 
directly proportional to a^(v)v, where o*(v) is the charge-exchange 
cross section. Hence this probability is 

o*(v)vf d9 d(P dv 

^2lT TSIt ^ 

j . I I o*(v)vf d6 d(j) dv 

M] M) M) 

where the triple integral is required for renormalization. Combining 
equations (Bl) and (B2) gives the probability as 


-Mv^ / 2kT , 


cr^(v)e 


i 3 

V sin 6 d0 d<|) dv 
-Mv2/2kTj^ 

o^(v)e v^ sin 6 d9 d(|i dv 


Equation (B3) is more tractable if the numerator is transformed from 
polar coordinates to cylindrical coordinates (v^.ttijVr), where both coordi- 
nate systems are shown in the following sketch. 
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Line of sight 



Making this transformation and integrating (B3) over all values of 
and (j) gives the probability of a charge -exchange neutral having v 
in dv as 

z 



(B4) 


Transforming the upper integral from and to v and v^ and 

evaluating the second integral in the denominator gives the probability 
for Vg ^ 0 as 



-Mv^/2kTj^ 

c. (v)e dv 


(B5) 


'0 
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We assume A emits a photon in a discrete transition. If we neglect 
the relativistic term in the Poppler relation (ref. 14) then 


/ 


A = 


X-v 
0 z 


(B6) 


where Xq is the wavelength of the line center, c is the velocity of 
:j.ight, and A is the difference between the wavelength and Xq. Con- 
verting equation (B5) from Vg to A and from v to energy E gives 
the probability 


where 



Y 


2 2 
Me A 


2X 


(B7) 


(B8) 


Since the line shape function is normalized to 1 (see eq. (1)), expres- 
sion (B7) is equal to F dA. Although equation (B7) was derived for pos- 
itive A, the same result is also obtained for negative A. 

A special case of equation (B7) is obtained by letting be a 

constant. This gives 

Y^^i A- rr\ 



where erfc is the complenjentary error fxjnction* This is not the 
Gaussian profile usually encountered in Doppler broadening. 

The second moments of the two line profiles can be obtained by sub- 
stituting the respective FU into equation (4). For the a^(E) case, 
equations (4) and (B7) give 




( 6 ) 


E dE 
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If is known, may be obtained from equation (6) by numerical 

methods, for the constant case equation (6) gives 




kT 

KiiAo 


Me 


( 5 ) 


which gives in closed form. 
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APPENDIX C 

DERIVATION OF SECOND MOMENT CORRECTION DUE TO MAGNETIC 
SPLITTING AND FINE STRUCTURE 

In order to determine ion temperature from equation (5) or (6) , two 
corrections should be made to the second moment of the charge- 

exchange neutral line observed by scanning with a monochromator. Because 
it is independent of the line shape, the first correction to make is the 
magnetic splitting and fine structure correction. A line may or may not 
have fine structure. It is well known that a magnetic field generally 
splits a fine structure component (or line if there is no fine structure) 
into several components . In mirror devices these components are usually 
sufficiently broadened that they overlap, especially for char ge -exchange . 
neutrals. This makes the observed line appear broader than it would 
otherwise appear and makes the observed second moment larger. The 

correction to derived below applies for any degree of component 

overlap, any type coupling in the atomic structure, and any strength mag- 
netic field. However, it is assumed that the components have the same 
shape and second moment, although they will, in general, have different 
integrated intensities (the integrated intensity is the spectral light 
intensity of a component integrated over all wavelengths) , We neglect 
polarization, although if the magnetic splitting correction ever were 
large compared to it would have to be taken into account (grating 

monochromators are notorious polarizers) . 

We assume each component has the same line shape function F nor- 
malized by 


r m - . 1 


(Cl) 


where is the wavelength of the center of the component, Also^ each 

of the k components has an integrated intensity normalized by 


E 

- 1=1 



(C2) 


where i -specifies ^ich component. The intensity of a component i at 
X is proportional to Ij^F(X - Xoi) • center of the k components 

collectively is 
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Changing the order of integration and summation gives 



(C3) 


(C4) 


But 


s 


XF(X - XQ.)dX 5 Xq^ 


(C5) 


Combining equations (C4) and (C5) gives 


K 


^Oi 


(C6) 


i=l 


The observed second moment is 







(X - Xo)“ ^ IiF(X - Xo^)dX 
i=l 



(C7) 


Y2 IiF(A - Xoi)dX 
i=l 


which can.be rewritten 


- ^Oi) + (^Oi - - Aoi)dX (C8) 


'0 



! 
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To simplify (C8) we define the second moment of component i by 




2 


X 



(A - AqP^ F(A - XQ^)dA 


(C9) 


Since it has the same value for all k components, no subscript i is 
needed on \^^)2 (the 2 does not refer to a component). Combining equa- 
tions (C8) and (C9) gives 


k 

i=l 

Solving for 

k 

i=l 


For any atomic line the I^^'s and the quantities Aq^ - Aq can be calcu- 
lated in principle if the magnetic field strength, direction of observa- 
tion relative to the magnetic field, and fine structure are known. For 
helium and hydrogen the only relevant cases are (1) Zeeman effect at low 
field strength, (2) intermediate coupling at intermediate field strength, 
and (3) Paschen-Back effect at high field strength. 


After the fine structure has been removed from the line shape (cor- 
responding to the use of eq. (Cll)), the residual line shape may be 
assumed symmetrical about Aq. 
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APPENDIX D 

DERIVATION OF SECOND MOMENT CORRECTION DUE TO 
SLIT FUNCTION OF MONOCHROMATOR 

After making the correction for magnetic splitting and fine struc- 
ture described in appendix C, the effect of the slit function of the 
monochromator (ref. 35) should be taken into account in the second moment 
When a line is scanned with a monochromator, the line appears wider than 
it really is due to the finite width of the entrance and exit slits, 
optical aberrations, and unwanted diffraction. These effects are ex- 
pressed collectively by the slit function of the monochromator, which is 
proportional to the response of the monochromator when set at one wave- 
length but illuminated by a different wavelength. Correcting for the 
slit function always reduces the second moment of a line shape. 

We assume the slit function S and the line shape function F are 
symmetrical (any fine structure having been removed by application of 
equation (Cll) to the second moment and a corresponding decomposition of 
the line shape function - although in practice the decomposition does not 
have to actually be carried out if only the temperature is desired) , Let 
be the wavelength at which the monochromator is set and A be the 
wavelength of some monochromatic light. We define 

X = A - A (Dl) 

m 

The slit function of many monochromators at medium to high resolution can 
be approximated by a Gaussian: 

S(x) =1 g-4x^ln2/w2 

w V “n* 

where W is the full half width at half intensity in wavelength units 
and S(x) is normalized by 


S(x)dx = 1 


The second moment 




of S(x) is 


0 - 




8 In 2 


(D3) 


(D4) 


Hence if W is measured experimentally, the second moment can be found 
from equation (D4) . The line shape function Fg observed at the exit 
slit as the spectrum is scanned is obtained from Lincke (ref. 35) after 
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considerable manipulation. 


m 




F(X - Xq)S(X - X^)dX 


(D5) 


where F(X) is the true line shape function of the line incident upon the 
entrance slit. Because S(X - X™) “ S(X™ - X) equation (D5) can be re- 
written 


F 

e 


m 



(D6) 


whe re 


A 


m 



(D7) 


which agrees with Griera (ref. 15). Mathematically Fq is simply the 
convolution of F and S. According to Papoulis (ref. 36), the charac- 
teristic function $g(o)) of Fg is given by 

$g(co) = ^p(aj)$g(oj) (D8) 


where $p(w) and $g(m) are the characteristic 
respectively. The second moment (or variance) 
(ref. 36) 


functions of 
(a 2)2 of Fg 


F and S , 
is then 



(D9) 


where in Identifying V^^/2 with Fg we have assumed that the line inci- 
dent on the entrance slit has had the splitting removed from it (appen- 
dix C) ■ 


The characteristic function of F can be found from equation (BIO) 
and is 


«p(«) “ 2 ■ 4A 2 ' ~ 4A 


) + J f= - £2) 


where is a degenerate hypergeometric function (ref. 37) and 
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A. -“si. 

2kTXp 



From equations (D2) and (D4) 

. . -l/2(x^)m^ 

$g(u)) = e 

From equations (D9) , (DIO), and (D13) and reference 37 



From equations (B12) , (D4) , (Dll) , (D12) , and (D14) 

<^^>3 - <^^2 - 

This is the same law that applies if both F and S are Gaussian. 


(Dll) 

(D12) 

(D13) 

(D14) 

(D15) 
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TABLE I. - SECOND MOMENT OF CHARGE-EXCHANGE NEUTRAL 
COMPONENT OF TWO LINES AS A FUNCTION 


OF ION TEMPERATURE 


Ion temperature^ 
Ti> 
eV 


Hydrogen Balmer 
alpha (Ha) 

(a) 

Helium 587 • 6 nm 
line 

20 

0.4429x10-1 

0.5817x10-2 

25 

.5228 

.6356 

30 

.5990 

,6885 

40 

.7437 

,7922 

50 

.8807 

.893 7 

60 

0.1012xloO 

0.9934 

70 

.1139 

.1092x10-1 

80 

.1262 

.1189 

100 

.1495 

.1382 

140 

.1914 

.1768 

200 

0.2450 

0.2352 

240 

.2766 

.2746 

300 

. 3206 

.3346 

400 

.3891 

.4361 

500 

.4562 

.5394 

i 

600 

0.5243 

0.6452 

700 

.5940 

.7547 

800 

.6652 

.8680 

1 .000 

,8106 

i ,1104x10® 

1 400 

.1104xl0l 

! .1580 

1 

2 000 

0.1539 

1 

0.2255 

2 400 

.1827 

.2689 

3 000 

.2260 

.3343 

4 000 

.2990 

,4482 

5.. 000 

.3718 

.5692 

6 000 

0.4424 

0.6954 

7 000 

.5093 

.8235 

8 000 

.5719 

.9499 

10 000 

.6837 

,1190xlol 

14 000 

.8631 

.1604 


^2 background gas assumed. 
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Figure 1. - Test section apd magnetic MeW configuration. 
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Figure Z - Eiecirode configuration and magnetic field lines. 



Figure 3. - Apparatus for Doppler broadening, shift, and line ratio measurements 
(not to scale). 
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Figure 4» - Comparison of various He 587,6 
nanometers line shapes for a He'*' tem- 
perature of 3000 electron volts. 


12x10^ 



ION TEMP ASSUMING ALL 
BACKGROUND, eV 

CS-7142B ^ 


Figure 5. ■■ Comparison of H'*' tem- 
peratures from Hj^assuming 
different background gases. 
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Figure 7. - Comparison of Ion temperatures 
from Hcc and Hp assuming H 2 background. 
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Figure 10. - H"*" bn temperature from Hr versus P/pB. 
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Figure 11. ^ Comparison of H'*' temperatures 
from Hr from neutral particle analyzer. 
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Figure 13. - Comparison of He'*' temperatures 
from He 587,6 and He 667. 8 nanometer lines. 
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Figure 14. - He'*' ion temperature 
from He 587.6 versus P/B. 
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Figure 15. - He"^ ion temperature from He 587. 6 
versus P/pB. 
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Figure 16. - Comparison of He“^ temper- 
atures from He 587.6 nanometer line 
and from neutral particle analyzer. 
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Figure 17, - Shift of wide 
component versus chord 
position for H«.. 
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Figure 18. - Electron temperature 
versus ratio of electrode voltage 
to background pressure for He. 
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